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. Abstract. We discuss the application of an effective field theory (EFT) which incorporates the 

chiral symmetry of QCD to Compton scattering from the proton and deuteron. We describe the 
chiral EFT analysis of the 7p scattering database presented in our recent review [1], which gives: 

ail = 10.5 ±0.5(stat) ±0.8(theory) , ^'^^ ^ 2.7 ±0.5 (stat) ±0.8 (theory), 
O , for the electric and magnetic dipole polarizability of the proton. We also summarize Ref. [l]'s chiral 
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EFT analysis of the world data on coherent Compton scattering from deuterium, which yields: 
ail = 10.5 ±2.0(stat) ±0.8 (theory) , j3|^'| = 3.6 ± l.O(stat) ±0.8(theory). 
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cn : THEORETICAL BACKGROUND 

O 

^ ■ Experiments to measure Compton scattering from the proton and deuteron are presently 

being pursued at a number of facilities around the world, including MAMI (Mainz) [2, 
3], HlyS at TUNL [4], and MAX-Lab at Lund [5]. Chiral effective field theory (;fEFT) is 
one of the main theoretical techniques used to analyze these experiments. xEVI gener- 
ates the most general Compton amplitude that is consistent with electromagnetic gauge 
5^ ! invariance, the pattern of chiral-symmetry breaking in QCD, and Lorentz covariance, to 

any given order in the small parameter P = {(U,m;r}/A, with co the photon energy, nin 
the pion mass and A the breakdown scale of the theory. 

The pioneering calculations of )p scattering in ;tEFT [6, 7] were performed in a 
theory with only nucleons and pions as explicit degrees of freedom. This reduces the 
breakdown scale A to the energy at which the A(1232) is excited, i.e. M/^ — ~ 300 
MeV. In Refs. [6, 7] nucleon scalar dipole polarizabiUties are predicted to be, at i^{P^), 

OCEI = 10/3mi = ^^^ = 12.6 X 10-4 fml (1) 



The corresponding cross sections agree well with experiment up to at least CO ~ nij,;, but 
do not capture the rise of the data towards the A(1232) peak. In this variant of ;tEFT 



even an ^(P^) calculation cannot describe the data at backward angles once O) > 180 
MeV [8, 9, 10]. 

The inclusion of the Delta as an active degree of freedom in the theory is therefore 
essential if the full power of the world's Compton data to shed light on fundamental 
hadron-structure parameters, e.g. polarizabilities, is to be realized. With the A(1232) 
included in ;tEFT [11, 12, 13] the ratio {M^ — M]\j)/A becomes one of the theory's 
expansion parameters. Ref. [14] pointed out that (Ma —Mn)/A is numerically rather 
similar to the expansion parameter in A-less calculations, nij^/ (Ma —M^), and denoted 
both as 5. In the "5-counting" adopted in Ref. [14] powers of the electronic charge e 
are shown explicitly, while Refs [6, 7] counted e ~ i°. Thus the Thomson amplitude 
is ^(p2) ff{e^b^) and structure effects start with ;cN loops at (^(P^) ~ ^(e^S^) in 
the low-energy region. Further, since Ma — M^ ~ 5, whereas m-n ~ 5^, %h. loops are 
suppressed by an additional power of 5, and do not enter the amplitude until G{e^b^^. 

The Delta-pole graph has a special role in 5-counting: it too is &{e^^^) for o ~ m^j:, 
but it becomes enhanced in the region © ~ Ma — Ma?, because of proximity to the Delta's 
on-shell point [14]. In this domain the effects that generate the resonance's finite width 
must be resummed. The dominant scattering mechanism in this "medium-energy" 
region is then the excitation of a dressed A(1232) by the magnetic transition from the 
nucleon state, followed by de-excitation via the same Ml transition. This effect occurs at 
At G{e^?>^) the B2N^ A(1232) transition must also be considered. Further 
discussion of ;tEFT, and this power counting, can be found in Ref. [1]. 

Eq. (1) is also the prediction for the neutron polarizabilities in (Delta-less) ;^EFT at 

0'(P^). To access a"^l and j8^j experimentally a nuclear target is required. The deuteron 
is the simplest nucleus, and an accurate description of its structure is obtained when 
XEPT is applied to the NN problem (see, e.g., Refs. [15, 16]). yd scattering was first 
calculated in (A-less) ;fEFT in Ref. [17] at &{P^), where it was demonstrated that there 
are large isoscalar yNN mechanisms ("exchange currents") which are crucial to 

obtaining reasonable agreement with the data. This calculation was extended to higher 
orders, and augmented with A(1232) degrees of freedom, in Ref. [9, 10, 18]. However, 
it was not until the work of Ref. [19] that a ;fEFT treatment of Compton scattering from 
deuterium which respected the Thomson limit for the yd amplitude was formulated. This 
computation also included A( 1232) degrees of freedom, and demonstrated the ability of 
;i^EFT to describe deuterium Compton data from threshold to ft) > 100 MeV. 

A NEW ANALYSIS OF yp SCATTERING IN ;tEFT 

The calculation of the yp differential cross section presented here includes the nucleon 
Born graph and the r-channel TU^ pole graph (both calculated covariantly). The A-pole 
graphs (s- and M-channel) are dealt with as described in Refs. [1, 14, 20] — covariantly 
and with a finite width stemming from ;rN loops. Compton ;rN and nA loop graphs are 
also added, so the amplitude includes effects which are of leading or next-to-leading or- 
der throughout the kinematic region < ftJiab ^ 350 MeV (apart from the loop correction 
to the yNA vertices), and all effects up to next-to-next-to-leading order — 6'{e^5^) — in 
the low-energy region ft) ~ m;,;. In addition, we include the contact interactions which 



encode the short-distance (r ^ l/nij^) contributions to the scalar polarizabilities. Strictly 
speaking these do not occur until 0{e^5^) in the low-energy domain, but they are nec- 
essary for an accurate description of yp data in HB;|^EFT (c.f. the case where pion loops 
are calculated relativistically [20, 21]). The coefficients of these contact terms are fit to 

the TP database, which is equivalent to fitting a^^ and pj^l . 

The parameters used in the ;rN sector take standard values (see Ref. [1]). The A( 1232) 
parameters — Mf^ — 293 MeV and gjiNf;, — 1 .425 are obtained from the Breit-Wigner 
peak and width, the latter via the relativistic formula. We adopt ^2/^1 = —0.34 for the 
ratio of E2 and Ml couplings (c.f. the ;|^EFT study of pion photoproduction [22]). 

Three EFT parameters — bi, Oi"^}, and — remain to be determined. They are fit to 
the yp data base discussed extensively in Ref. [1]. For the reasons explained there, in 
the low-energy region we include data from Refs. [23, 24, 25, 26, 27, 28, 29, 30, 31]. 
We float the normalization of each of these data sets within the quoted normalization 
uncertainty. In the medium-energy region the data sets of Refs. [29, 32] are in significant 
disagreement with those from MAMI (most notably Refs. [33, 34]) and a consistent fit 
cannot be obtained. We have chosen to use the MAMI data for our fits in this region [1]. 

Our strategy is to determine the 7NA Ml coupling b\ by considering the MAMI data 

for (0[ab =200-325 MeV, then fit a^^} and j8^} to the low-energy data up to 170 MeV, 
and iterate until convergence is reached. The j^/d.o.f. of the low-energy Hallin data 
is hard to accept, and so we prefer to quote our best results without them. We then 
obtain a solution with a ;t:^/d.o.f. = 106.1/124. This "modified ff{e^8'^)" fit yields 

b\ = 3.66 ±0.03 and the values of and jij^l quoted in the abstract. Since these sum 

to a value consistent with the Baldin sum-rule constraint a'^^ + pj^l = 13.8 ±0.4 [31], 

one can impose this relation to find a^^ - ^j^l = 7.7 ±0.6 with ;t;^/d.o.f. = 106.5/125, 
an unchanged bi, and the cross sections displayed in Fig. 1. These fits are stable against 
reasonable variations in the procedure [1], and agree with data [33, 34, 35, 36, 37, 38, 39] 
well beyond the region in which the free parameters are determined (see Fig. 1). 



A NEW ANALYSIS OF yd SCATTERING IN ;tEFT 

The ;|^EFT treatment of yd scattering developed in Ref. [19] is valid from threshold 

to cOiab ^ 100 MeV, and represents a complete (modified) ^{e^5^) calculation in both 
the two- and one-nucleon sectors. It has the added virtue that the dependence of cross 
sections on the choice of deuteron wave function is < 1% (c.f. Ref. [10]). We now 
employ the ;tEFT deuteron wave function at NNLO [15] within this theory to extract 

a^l and jS^^'j from the yd elastic data of Refs. [40, 41, 42]. 

The fit to these isoscalar combinations of scalar dipole polarizabilities yields the 
results quoted in the abstract, with ;^^/d.o.f. = 24.3/24 (see Fig. 2). In contrast to the 
proton case, this is a consistent data base: each experiment contributes roughly equally 
to the Slid the extracted polarizabilities are largely insensitive to the elimination of 
any one data set. The isoscalar polarizabilities we obtain are close to the proton ones, so 
isovector effects in agi and /3mi are small, as predicted by ;tEFT at ff{P^). We also used 
the (isoscalar) Baldin constraint to reduce statistical uncertainties in a one-parameter fit. 
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FIGURE 1. Comparison of our ;i;EFT result for yp scattering with data. The lab cross section in nb/sr is 
shown in 10° bins for 0iab as a function of (Oiab in MeV, with insets showing the fit region. The grey band 
shows the variation within the statistical error of the one-parameter fit. Adapted from Ref. [1], where the 
legend for experimental data can be found. 
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FIGURE 2. yd cross sections at 49 and 94.5 MeV in the two-parameter (dashed) and one-parameter 
(solid) determinations of the isoscalar spin-independent dipole polarisabilities. Bands: statistical error of 
the Baldin constrained fit. Data at 49 (94.5) MeV from Ref. [40] ([41]). Adapted from Ref. [1]. 



and obtained very similar results (see Fig. 2). For further details, see Ref. [1]. 
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